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Abstract: Pentatomomorpha is one of the most biodiverse infraorders among the true bugs (Hemiptera: 
Heteroptera). Phylogenetic relationships among the superfamilies within this infraorder have been uncertain, 
especially for the Eutrichophora. The previous studies were based on morphological characters, or just 
mitochondrial or nuclear genes, or only partial 18S rDNA and COL. In this study, we used maximum likelihood 
(ML) and Bayesian inference (BI) based on massive molecular datasets (18S rDNA, 28S rDNA, Hox and 
mitochondrial genes totaling 21 loci and 12,538 characters) to infer a robust phylogeny for this terrestrial group. 
Results strongly support the monophyly of all superfamilies; the superfamily status of Aradoidea and the 
following relationships: (Aradoidea + (Pentatomoidea + (Coreoidea + (Lygaeoidea + Pyrrhocoroidea)))) in 
Pentatomomorpha, and (Coreoidea + (Lygaeoidea + Pyrrhocoroidea)) in Eutrichophora. Our results suggest that 
sampling greater numbers of genes is an effective tool for resolving phylogenetic problems. 
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Introduction 


Pentatomomorpha (Hemiptera: Heteroptera) includes some of the largest terrestrial bugs. 
This group causes economic damage to many species of crops because its members are chiefly 
phytophagous (Schuh & Slater 1995). Pentatomomorphan bugs comprise more than 14,000 
species (Weirauch & Schuh 2011). This large number of species has been variously grouped 
into four (Schaefer 1964), five (Stys 1961; Schaefer 1993; Schuh & Slater 1995), six (Henry 
1997), or seven (Schuh 1986; Henry & Froeschner 1988) superfamilies (Table 1) and 39 
families (Schuh & Slater 1995; Henry 1997; Rider 2015). The superfamily Aradoidea was 
raised to the infraorder Aradomorpha by Sweet (1996, 2006) based on morphological evidence, 
together with ecological, fossil and biogeographical evidence, but this viewpoint was treated 
as invalid in Cassis and Schuh (2010). 


Table 1. Previous classification systems of Pentatomomorpha 


Schaefer 1964 Stys 1961 Schaefer 1993 Schuh1986 Henry 1997 
Schuh & Slater Henry & 
1995 Froeschner 1988 

Aradoidea Aradoidea Aradoidea Aradoidea Aradoidea 

Pentatomoidea Pentatomoidea Pentatomoidea Pentatomoidea Pentatomoidea 

Coreoidea Coreoidea Coreoidea Coreoidea Coreoidea 
Pyrrhocoroidea Pyrrhocoroidea Pyrrhocoroidea 
Lygaeoidea Lygaeoidea Lygaeoidea 

Piesmatoidea Piesmatoidea Piesmatoidea 

Idiostoloidea Idiostoloidea Idiostoloidea 


Phylogenetic relationships proposed previously at the superfamily level within 
Pentatomomorpha are summarized in Fig. 1. The works of Henry (1997) (Fig. 1A) and Yao et 
al. (2012) (Fig. 1E) were based on morphological characters. Li et al. (2005) (Fig. 1B) 
analyzed the phylogeny of Pentatomomorpha using molecular sequence data (partial 18S 
rDNA and COI mtDNA) for the first time, and supported the monophyly of Pentatomomorpha, 
a sister-group relationship of Aradoidea and Trichophora, and a clade comprising Lygaeoidea, 
Coreoidea, and Pyrrhocoroidea. Later, several other works based on molecular sequences were 
published (Hua et al. 2008; Tian et al. 2011; Yuan et al. 2015). Of these, Hua et al. (2008 ) 
(Fig. 1C) and Yuan et al. (2015) (Fig. 1D) were both mitochondrial genome sequences-based 
analyses, but reached different relationships for the Eutrichophora, a concept proposed in Xie 
et al. (2005). In Tian et al. (2011) (Fig. 1D), 26 pentatomomorphan and four cimicomorphan 
putative families and six Hox gene fragments were analyzed, and the relationship (Aradoidea 
t (Pentatomoidea + (Lygaeoidea + (Coreoidea + Pyrrhocoroidea)))) was supported. 

Until now, the monophyly of Pentatomomorpha, and the close relationship between 
Pentatomoidea and Eutrichophora have been generally accepted (Henry 1997; Xie et al. 2005; 
Li et al. 2005; Hua et al. 2008; Tian et al. 2011; Yao et al. 2012; Yuan et al. 2015). The 
phylogenetic relationships within the Pentatomomorpha, however, have not reached full 
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agreement. The unsolved problems are: 1) the monophyly of the superfamilies (Xie et al. 2005; 
Li et al. 2005); 2) whether Aradoidea should be treated as superfamily Aradoidea or infraorder 
Aradomorpha (Sweet 1996, 2006); and 3) the phylogenetic relationships among the 
superfamilies, especially in the Eutrichophora (Fig. 1). 

Previous phylogenetic studies were mostly based on morphological characters, fossils or 
several single source genes. Lack of a sufficient number of features, errors in accessing 
homology among features of various organisms, and convergent evolutionary changes can 
affect the accuracy of phylogenetic inference (Boore & Brown 1999). Multiple genes with 
increased sequence lengths are favorable for accurate phylogenetic analyses (Wortley ef al. 
2005; Kawahara et al. 2011; Corl & Ellegren 2013). Sequences used to construct phylogenetic 
trees of Pentatomomorpha were less than 1.5 kb (Li et al. 2005), or just from nuclear (Tian et 
al. 2011) or mitochondria (Hua et al. 2008; Yuan et al. 2015). In this study, we took advantage 
of nucleotide and genomic sequencing projects to assemble a data set of 21 genes and 37 
species. Based on this data set, a preliminary phylogenetic framework of Pentatomomorpha is 
proposed and the status of Aradoidea is analyzed. 


A B 
Aradoidea 


Pentatomoidea 
Aradoidea 


Pyrrhocoroidea 
Pentatomoidea 


Coreoidea Lygaeoidea 
+ 
Idiostoloidea Pyrrhocoroidea 
+ 
Lygaeoidea Coreoidea 
C D E 
Aradoidea Aradoidea 
. Aradoidea 
Pentatomoidea Pentatomoidea 
Pyrrhocoroidea Lygaeoidea Pentatomoidea 
Lygaeoidea Pyrrhocoroidea 
Coreoidea Coreoidea opa 


(sensu lato) 


Figure 1. Proposed phylogenetic hypotheses within the Pentatomomorpha. A. Henry (1997); B. Li et al. 
(2005); C. Hua et al. (2008); D. Tian et al. (2011), Yuan et al. (2015); E. Yao et al. (2012). 


Material and methods 


Taxa sampling. Taxa were chosen to represent the major lineages according to the 
classification of Schuh and Slater (1995), Henry (1997) and Rider (2015). The terminal taxa 
were decided according to Campbell and Lapointe (2009). In total, 37 species were sampled, 
of which 28 species were ingroups. Two gerromorphans, two nepomorphans, two 
leptopodomorphans and three cimicomorphans were selected as successive outgroups (Table 2). 

Gene sampling. A total of 15 fragments of Hox genes of Ptilomera tigrina Uhlen, 1860, 
Microvelia douglasi Scott, 1874, Tingis ampliata (Herrich-Schaeffer, 1838) and Anthocoris 
thibetanus Poppius, 1909 were initially sequenced in this study. All of the other Hox genes, all 
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18S rDNA and 28S rDNA sequences, and most mitochondrial genome genes were generated 
by authors of the previously published studies (Hua et al. 2008; Tian et al. 2011; Li et al. 2012; 
Wang et al. 2015) (Table 2). 

Laboratory procedure and primers. The laboratory procedures were the same as those 
described in Li et al. (2012). Our newly amplified five Hox gene exons are: the 5' end of 
abd-A, 3’end of Dfd, 5’ end of Dfd, 5’ end of Ubx and 3’ end of pb, with sequence lengths 
ranging from 450 to 1,000 bp. Universal degenerate primers were used in this study. Some are 
from Tian et al. (2011) and Li et al. (2012); some were newly designed as shown in Table 3. 


Table 3. Universal degenerate primers used in this study 


Upper primer (5’-3) Lower primer (5’-3) 
5’ end of abd-A Tian et al. 2011 
3' end of Dfd Tian et al. 2011 
5' end of Dfd Tian et al. 2011 
5' end of Ubx Li et al. 2012 
3' end of pb Tian et al. 2011 
3' end of pb for GTRAA RGTYT GGTTY CAAAA YCG TAGTA YTCSG GCGYR AAGTC GTT 
Gerromorpha 
5' end of abd-A TACCANCCNCAAATGCTGCC TGTCWGGCATWGGTTGKGAAG 


for Gerromorpha 


Sequences of 18S rDNA, 28S rDNA, Hox and 13 protein coding genes of mitochondrial 
genomes were aligned by Muscle, embedded within MEGA 6.06 (Tamura et al. 2013) under 
the default conditions of the program. The alignments of 18S rDNA and 28S rDNA were 
adjusted according to the secondary structure models of heteropteran 18S rRNA and 28S 
rRNA (Xie et al. 2009; Wang et al. 2013), and length-variable regions of both genes were 
eliminated before phylogenetic reconstructing (Wang et al. 2015). The six Hox genes were 
also manually corrected according Tian et al. (2011). Finally, 18S rDNA (1937 bp), 28S rDNA 
(4163 bp), six Hox genes (2708 aa) and 13 mitochondrial protein coding genes (PCGs) (3730 
aa) were selected as molecular markers. All of the four data sets (18S rDNA, 28S rDNA, Hox 
and mt-genomes) were concatenated in DAMBE V5.3.48 (Xia 2013). 

Phylogenetic analyses were conducted utilizing maximum likelihood (ML) and Bayesian 
inference (BI). Best-fitting substitution models were selected by Modeltest 3.7 (Posada & 
Crandall 1998) and Treefinder (Jobb 2011) for the nrDNAs and amino acids, respectively. ML 
analyses were conducted using RAxML 8.0.12 in PThreads version (Stamatakis 2014). The 
most appropriate substitution model was GTR+G+I for 18S rDNA and 28S rDNA, mtArt + G 
+ I for ATP6, COI, COII, COII and ND4L, mtZOA + G + I for ATP8, CYTB, ND1, ND2, 
ND3, ND4, NDS, ND6 and 5’ end of Scr, Dayhoff + G + I for 5’ end of abd-A and 3’ end of 
Dfd, JTT + G+ I for 5’ end of Dfd and 5’ end of Ubx, MIX [Dayhoff, BLOSUM, JTT, LG 
PMB, VT, WAG, mtREV, mtMam, mtArt, mtZOA, cpREV, rtREV, betHIV, witHIV, FLU] + G 
+ I for 3’ end of pb. The best ML tree was calculated by RAxML with 100 runs, followed by 
2000 bootstrap replicates. The consensus tree was built based on the best ML tree and the 
bootstrap tree. 

Bayesian analyses were conducted using MrBayes 3.2.5 (Ronquist et al. 2012) on CPUs 
with the gaps in the matrices treated as missing. The GTR substitution model for 188 rDNA 


Entomotaxonomia (2016) 38(2): 81-93 85 


and 28S rDNA, and the mixed substitution model for the amino acids were used under a 
gamma distribution (+G) with a proportion of invariable sites (+I) to account for among-site 
rate variation. Four chains were used and the datasets were run for 5,000,000 generations until 
the average standard deviation of split frequencies between the two-parallel Markov chain 
Monte Carlo runs dropped below 0.01. Trees were sampled every 1000th cycle, with a burn-in 
of 25% of sampled trees. 


Phylogenetic results 

Two trees generated through two methods (ML and BI) recovered an identical topology 
for the higher-level nodes (Figs. 2, 3). The monophyly of all superfamilies except Aradoidea of 
Pentatomomorpha was supported. Aradoidea was the sister group of the remaining 
pentatomomorphan taxa, and then Pentatomoidea came out. The relationships among 
Eutrichophora were as follows: (Coreoidea + (Lygaeoidea + Pyrrhocoroidea)). The suggested 
relationships among the nine outgroups and Pentatomomorpha in this paper were (Gerridae + 
Veliidae + (Nepomorpha + (Leptopodomorpha + (Cimicomorpha + Pentatomomorpha)))). 
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Figure 2. ML phylogram inferred from the combined data. Bootstrap support values are indicated at each node 
(> 50%). 


Two trees generated through two methods (ML and BI) recovered an identical topology 
for the higher-level nodes (Figs. 2, 3). The monophyly of all superfamilies except Aradoidea of 
Pentatomomorpha was supported. Aradoidea was the sister group of the remaining 
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pentatomomorphan taxa, and then Pentatomoidea came out. The relationships among 
Eutrichophora were as follows: (Coreoidea + (Lygaeoidea + Pyrrhocoroidea)). The suggested 
relationships among the nine outgroups and Pentatomomorpha in this paper were (Gerridae + 
Veliidae + (Nepomorpha + (Leptopodomorpha + (Cimicomorpha + Pentatomomorpha)))). 

ML tree (Fig. 2). Branch support for the monophyly of all the superfamilies in 
Trichophora was robust (BP = 100%), as well as for the basal branch, and clade (Coreoidea + 
(Lygaeoidea + Pyrrhocoroidea)). The support indices of clades (Pentatomoidea + (Coreoidea + 
(Lygaeoidea + Pyrrhocoroidea))) and (Lygaeoidea + Pyrrhocoroidea) were 91% and 55%, 
respectively. 
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Figure 3. BI phylogram inferred from the combined data. Bayesian posterior probabilities values are indicated 
at each node (> 50%). 


BI tree (Fig. 3). Bayesian topologies were identical to those recovered in the ML tree (Fig. 
2). The posterior probability values supporting the monophyly of all superfamilies and the 
relationships among the superfamilies were all 10096. In addition, inter-family relationships 
were also well-resolved. 


Discussion 


Pentatomomorpha is recognized as a monophyletic group by most researchers. Our 
results were in accordance with this opinion with high support values in ML and Bayesian 
trees. 

The monophyly of Pentatomoidea and the sister relationship of Pentatomoidea and the 
other Trichophora were strongly supported in both analyses (BP > 90%; PP = 1.00), views that 
corroborate almost all of the recent studies (Henry 1997; Li et al. 2005; Hua et al. 2008; Yao et 
al. 2012; Yuan et al., 2015). Although only one representative was chosen for Aradoidea, no 
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one has questioned its monophyly up to now. So the monophyly of Aradoidea is also accepted 
in this research. 

The superfamily Aradoidea was raised to the infraorder Aradimorpha by Sweet (1996, 
2006) based on morphological, ecological, fossil and biogeographical evidences. In the 
mitochondrial genome analyses of Hua et al. (2008), they found a unique gene rearrangement 
in Aradoidea as its autapomorphy. However, Aradoidea was consistently placed as an 
independent branch of Pentatomomopha in our research. The sister relationship of Aradoidea 
and Trichophora (Leston 1958; Henry 1997; Li et al. 2005; Hua et al. 2008; Yao et al. 2012; 
Yuan et al. 2015) was further confirmed in this research. 

As for the monophyly and relationships of the other three superfamilies i.e. Coreoidea, 
Lygaeoidea and Pyrrhocoroidea, there were different views: Henry (1997) gave the topology 
((Coreoidea + Pyrrhocoroidea) + Lygaeoidea) based on morphological characters. Li et al. 
(2005) recognized them as a paraphyletic group based on partial 188 rDNA and COI 
sequences; Xie et al. (2005) performed a Bayesian analysis using 18S rDNA and reported the 
hypothesis of (Pyrrhocoroidea + (Coreoidea + Lygaeoidea)); Hua et al. (2008) used mitochondrial 
genomes as molecular markers and obtained the same results as Xie et al. (2005). Tian et al. 
(2011) analyzed six Hox genes and found the same relationship as Henry (1997). Yao et al. 
(2012) clarified that the remainder of the Trichophora (except for Pentatomoidea) should be 
assigned to 5 superfamilies (Coreoidea, Pyrrhocoroidea, Idiostoloidea, Lygaeoidea, and 
Piesmatoidea) based on a combination of fossil and extant morphological characters, but 
didn't give the relationships among them. Phylogenetic analyses based on mitogenomic data 
of Yuan et al. (2015) supported a phylogenetic relationship of (Lygaeoidea + (Pyrrhocoroidea 
+ Coreoidea)), the same result found by Henry (1997) and Tian et al. (2011). 

In this study, the relationship ((Coreoidea + (Lygaeoidea + Pyrrhocoroidea)) was 
demonstrated. The three superfamilies are united by four synapomorphies with high 
consistency indices: presence of an m-chromosome, typically (ground plan) three trichobothria 
on at least some of the abdominal segments, 3-lobed (or more) salivary glands, and loss of 
dorsal abdominal scent gland in nymphs between tergites 3/4. Further, several studies have 
suggested that the incomplete abdominal sutures found in some taxa (e.g., Physopeltinae 
[Largidae] and some Pyrrhocoridae) link these bugs with the Rhyparochrominae in the 
Lygaeoidea (Henry 1997). 

For the relationship among outgroups, (Nepomorpha + (Leptopodomorpha + 
(Cimicomorpha + Pentatomomorpha))) was supported and this is in congruence with results of 
previous analyses of morphological data (Schuh 1979) and molecular data (Wang et al. 2015). 

In this paper, we focused primarily on the higher taxa of Pentatomomorpha. Inclusion of 
wider sampling to establish the relationships among lower taxa would be the expected next 
step. 
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Table 2 Taxa and Genes Used in this Analysis 


Infraorder-level 


Gerromorpha 
Gerromorpha 
Nepomorpha 
Nepomorpha 
Leptopodomorpha 
Leptopodomorpha 
Cimicomorpha 
Cimicomorpha 
Cimicomorpha 
Pentatomomorpha 
Pentatomomorpha 
Pentatomomorpha 
Pentatomomorpha 
Pentatomomorpha 
Pentatomomorpha 
Pentatomomorpha 
Pentatomomorpha 
Pentatomomorpha 
Pentatomomorpha 
Pentatomomorpha 


Pentatomomorpha 


Pentatomomorpha 
Pentatomomorpha 
Pentatomomorpha 
Pentatomomorpha 
Pentatomomorpha 
Pentatomomorpha 
Pentatomomorpha 
Pentatomomorpha 
Pentatomomorpha 
Pentatomomorpha 
Pentatomomorpha 
Pentatomomorpha 
Pentatomomorpha 


Pentatomomorpha 


Pentatomomorpha 


Pentatomomorpha 


Family-level 


Gerroidea:Gerridae 
Gerroidea: Veliidae 
Nepoidea:Nepidae 
Ochteroidea:Ochteridae 
Saldoidea:Saldidae 
Leptopodoidea:Leptopodidae 
Naboidea:Nabidae 
Miroidea:Tingidae 
Cimicoidea:Anthocoridae 
Aradoidea:Aradidae 
Pyrrhocoroidea:Largidae 
Pyrrhocoroidea:Pyrrhocoridae 
Pyrrhocoroidea:Pyrrhocoridae 
Lygaeoidea:Lygaeidae 
Lygaeoidea:Lygaeidae 
Lygaeoidea:Artheneidae 
Lygaeoidea:Blissidae 
Lygaeoidea:Cymidae 
Lygaeoidea:Geocoridae 
Lygaeoidea:Piesmatidae 


Lygaeoidea:Rhyparochromidae 


Coreoidea:Rhopalidae 
Coreoidea:Alydidae 
Coreoidea:Alydidae 
Coreoidea:Coreidae 
Coreoidea:Stenocephalidae 
Pentatomoidea:Tessaratomidae 
Pentatomoidea:Pentatomidae 


Pentatomoidea:Pentatomidae 


Pentatomoidea:Acanthosomatidae 


Pentatomoidea:Scutelleridae 
Pentatomoidea:Cydnidae 
Pentatomoidea:Plataspidae 
Pentatomoidea:Plataspidae 


Pentatomoidea:Urostylididae 


Pentatomoidea:Urostylididae 


Pentatomoidea:Urostylididae 


Species 


Ptilomera tigrina 
Microvelia douglasi 
Ranatra chinensis 
Ochterus marginatus 
Saldula saltatoria 
Valleriola buenoi 
Alloeorhynchus sp. 
Corythucha ciliata 
Anthocoris pilosus 
Mezira sp. 

Physopelta gutta gutta 
Dysdercus cingulatus 
Macrocheraia grandis 
Arocatus melanocephalus 
Sinorsillus piliferus 
Artheneis intricata 
Cavelerius excavatus 
Cymus koreanus 
Geocoris ater 

Piesma maculatum 


Neolethaeus assamensis 


Stictopleurus punctatonervosus 
Megalotomus ornaticeps 
Riptortus pedestris 

Cletus punctiger 
Dicranocephalus alticolus 
Eurostus validus 
Rhaphigaster nebulosa 
Eurydema maracandica 
Elasmucha laeviventris 
Cantao ocellatus 
Fromundus pygmaeus 
Coptosoma bifarium 
Paracopta maculata 


Urochela luteovaria 


Tessaromerus licenti 


Urostylis chinai 


Accession 
number of 
18SrDNA 
KJ461233 
KJ461281 
KJ461186 
KJ461251 
KJ461170 
KJ461310 
KJ461245 
KJ461201 
KJ461226 
KJ461285 
KJ461164 
KJ461263 
KJ461199 
KJ461273 
KJ461269 
KJ461299 
KJ461294 
KJ461232 
KJ461218 
KJ461171 
KJ461280 


KJ461217 
KJ461221 
AB725684 
KJ461173 
KJ461228 
KJ461181 
X89495 
JX997807 
KJ461262 
KJ461182 
KJ461296 
KJ461259 
KJ535895 
KJ461205 


KJ461291 
KJ461275 


Accession 


number of 


28SrDNA 
JF719825 

KJ461178 
JF719824 

KJ461315 
KJ461213 
KJ461290 
KJ461293 
KJ461175 
KJ461268 
KJ461189 
KJ461255 
KJ461235 
KJ461172 
KJ461308 
KJ461194 
KJ461203 
KJ461237 
KJ461312 
KJ461309 
KJ461169 
KJ461225 


KJ461286 
KJ461236 
AB725684 
KJ461219 
KJ461267 
JF719828 
EU426880 
JX997806 
KJ461208 
KJ461230 
KJ461209 
KJ461239 
KJ535890 
KJ461306 


KJ461300 
KJ461287 


Species 


Ptilomera tigrina 
Entomovelia sp. 
Laccotrephes robustus 
Ochterus marginatus 
Saldula arsenjevi 
Leptopus sp. 
Alloeorhynchus bakeri 
Corythucha ciliata 
Orius niger 
Neuroctenus parus 
Physopelta gutta 
Dysdercus cingulatus 
/ 


Kleidocerys resedae 


/ 
/ 
/ 


Geocoris pallidipennis 


" 


/ 

Stictopleurus subviridis 
f 

Riptortus pedestris 
Hydaropsis longirostris 
i 

Eusthenes cupreus 
Halyomorpha halys 
Eurydema gebleri 

/ 

/ 

Macroscytus gibbulus 
Coptosoma bifaria 


Megacopta cribraria 


Urochela quadrinotata 


/ 


Continued Table 2 on the right 


Accession 


number of Mt 


genomes 
KP400583 
KP400582 
NC_012817 
NC_012820 
NC_012463 
FJ456946 
NC_016432 
KC756280 
EU427341 
NC 012459 
NC 012432 
NC 012421 
/ 

KJ584365 


NC 012424 
/ 

/ 

NC 012888 
/ 

NC 012462 
NC 012456 
/ 

NC 022449 
NC 013272 
KP207595 

/ 

/ 

NC 012457 
NC 012449 
NC 015342 
JQ743678 

/ 

/ 


Species 


Ptilomera tigrina 
Microvelia douglasi 
Ranatra chinensis 
Ochterus marginatus 
Saldula pallipes 
Valleriola sp. 

Gorpis annulatus 

Tingis ampliata 
Anthocoris thibetanus 
Carventus hainanensis 
Physopelta quadriguttata 
Pyrrhopeplus posthumus 
/ 

Lygaeus equestris 
Oncopeltus fasciatus 
Artheneis intricata 
Dimorphopterus spinolae 
Cymus koreanus 
Geocoris pallidipennis 
Piesma josifovi 
Dieuches femoralis 
Aeschyntelus chinensis 

/ 

Riptortus pedestris 
Acanthocoris scaber 

/ 

Eurostus validus 

Nezara viridula 

y. 

Acanthosoma crassicauda 
Poecilocoris lewisi 
Aethus indicus 

/ 

Megacopta cribraria 
Urochela distincta 

/ 

/ 


Accession 
number of 
S'end abd-A 
0 

0 
JQ409392 
JQ409394 
JQ409389 
JQ409390 
FJ851727 
0 

0 
FJ851698 
FJ851712 
FJ851711 
y 
FJ851700 
FJ851728 
FJ851705 
FJ851703 
FJ851701 
FJ851704 
FJ851710 
FJ851707 
FJ851715 
/ 
FJ851713 
FJ851714 
/ 
FJ851722 
FJ851723 
/ 
FJ851717 
FJ851720 
FJ851718 
/ 
FJ851719 
FJ851716 
/ 

/ 


3'end Dfd 


0 

f 
JQ409405 
JQ409407 
JQ409402 
JQ409403 
FJ851754 
0 

0 
FJ851729 
FJ851739 
FJ851738 
/ 
FJ851731 
FJ851755 
FJ851734 
FJ851733 
/ 

/ 
FJ851737 
/ 
FJ851742 
/ 
FJ851740 
FJ851741 
/ 
FJ851749 
FJ851750 
/ 
FJ851744 
FJ851747 
FJ851745 
/ 
FJ851746 
FJ851743 
71 

/ 
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5'end Dfd 


0 

0 
JQ409432 
JQ409434 
JQ409429 
JQ409430 
FJ851784 
0 

0 
FJ851756 
FJ851769 
FJ851768 
f 
FJ851758 
/ 
FJ851762 
FJ851760 
/ 
FJ851761 
FJ851767 
FJ851764 
FJ851772 
/ 
FJ851770 
FJ851771 
/ 
FJ851779 
FJ851780 
/ 
FJ851774 
FJ851777 
FJ851775 
ff 
FJ851776 
FJ851773 
/ 

/ 


3'end pb 


0 

0 
JQ409445 
JQ409447 
JQ409442 
JQ409443 
FJ851810 
0 

/ 
FJ851785 
FJ851798 
FJ851797 
/ 
FJ851787 
/ 
FJ851792 
FJ851790 
FJ851788 
FJ851791 
/ 
FJ851794 
FJ851801 
/ 
FJ851799 
FJ851800 
/ 
FJ851808 
FJ851809 
/ 
FJ851803 
FJ851806 
FJ851804 
/ 
FJ851805 
FJ851802 
/ 

/ 


5'end Scr 


NN 7. 7. 7. CX 


FJ851840 
/ 

FJ851839 
FJ851811 
FJ851825 
FJ851824 
/ 

FJ851813 
FJ851841 
FJ851818 
FJ851816 
FJ851814 
FJ851817 
FJ851823 
FJ851820 
FJ851828 
/ 

FJ851826 
FJ851827 
/ 

FJ851835 
FJ851836 
/ 

FJ851830 
FJ851833 
FJ851831 
/ 

FJ851832 
FJ851829 
/ 

/ 


5'end Ubx 


/ 

/ 
JQ409419 
JQ409421 
JQ409416 
JQ409417 
FJ851870 
0 
FJ851868 
FJ851842 
FJ851856 
FJ851855 
/ 
FJ851844 
FJ851871 
FJ851849 
FJ851847 
FJ851845 
FJ851848 
FJ851854 
FJ851851 
FJ851859 
/ 
FJ851857 
FJ851858 
/ 
FJ851866 
FJ851867 
/ 
FJ851861 
FJ851864 
FJ851862 
/ 
FJ851863 
FJ851860 
/ 

/ 


Note: Sequences marked with 0 were newly amplified in this research and the accession numbers will 


be supplied after assigned by the GenBank staff. 


